Phenomenological Theory of Isotropic-Genesis Nematic Elastomers 
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We consider the impact of the elastomer network on the structure and fluctuations in the isotropic- 
genesis nematic elastomer, via a phenomenological model that underscores the role of network 
compliance. The model contains a network-mediated nonlocal interaction as well as a new kind 
of random field, which reflects the memory of the nematic order present at cross-linking, and also 
encodes local anisotropy due to localized polymers. Thus, we predict a regime of short-ranged 
oscillatory spatial correlations (both thermal and glassy) in the nematic alignment trapped into the 
network. 
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Consider a melt or solution of nematogenic polymers, 
by which we mean long, flexible polymers carrying rod- 
like units. These units, which give the system the pos- 
sibility of exhibiting liquid crystallinity, may be inte- 
grated along the polymer chain backbones (the main- 
chain case) or in groups that dangle from the backbone 
(the side-chain or pendant case). Now consider the pro- 
cess of instantaneous cross-linking. Here, one begins with 
the melt or solution at equilibrium and — so rapidly that 
hardly any relaxation has time to occur — one introduces 
permanent bonds between some random fraction of the 
pairs of chain segments that happen, at the instant of 
cross-linking, to be nearby one another. When the cross- 
linking process is carried out in the isotropic state of the 
nematogens, the resulting material is called an isotropic- 
qenesis nematic elastomer (or IGNE; see, e.g., Ref. [l|- 
|7[). It is a macroscopic random network medium that 
"memorizes" both the positions of the chain segments 
and the orientations of the nematogen units at the in- 
stant of cross-linking. The memorization is, however, 
only partial, as a result of the thermal fluctuations that 
occur in the new, post-cross-linking equilibrium state. 

To describe the structure and correlations of the sys- 
tem post cross- linking, we employ the local nematic order 
parameter Qdd'(Y), which is traceless, symmetric, and of 
rank- two, and is defined microscopically via 

p 

where P is the number of rod-like units, N'' is the mi- 
croscopic unit orientation vector of unit p and R^ is its 
microscopic position vector in £)-dimcnsional space, and 
r is an arbitrary position vector. In addition, we charac- 
terize the random local environmental anisotropy, which 
tends to induce local nematic alignment Q in the post 
cross-linking system, via the tensor Y^: 



Y,(r) = Y(r) + ^ j d^'r' H(r - r') Q"(r'). (2 



This random environmental anisotropy is caused by the 
thermally averaged part of random local spatial arrange- 
ment of the localized polymers at post-cross- linking equi- 
librium. It consists of two parts: (i) a part that is inde- 
pendent of the the pattern of local nematic alignment Q'^ 
present at the instant of cross-linking, which we call the 
memory-independent random field and denote by Y; and 
(ii) another part that is due to the pattern of Q°, which 
we call the memory- dependent random field. The pat- 
tern Q'' is partially imprinted in the network structure, 
and this imprint then partially elicits a response simi- 
lar to Q° in the post-cross-linking state. The relation- 
ship between Q and Q*' is characterized by a "smearing" 
kernel, which embodies the idea that Q (i.e., the post- 
cross-linking equilibrium-state memory of Q'^) is partially 
erased, as a result of the positional thermal fluctuations 
of the network. Equivalently, viewed from wave-vector 
space, the contribution from Q° becomes H,^Qk- Phys- 
ically, wc expect II(r) to be positive and bell-shaped, 
operative primarily over a region of order the typical lo- 
calization length ^i (which reflects how weakly localized 
the network constituents are; see, e.g., [8|), and to decay 
monotonically with increasing \r\ over this Icngthscale, 
ultimately tending to zero for \r\ 3> £,l- Correspond- 
ingly, in wave-vector space Hk would decay monotoni- 
cally to zero over a scale S,^^- Hence, we see that H 
serves as a "soft filter" de-amplifying — more strongly the 
shorter the lengthscale — the contributions made by the 
Fourier components of Q° to the random anisotropic en- 
vironment on distance scales shorter than S,l- This is a 
natural consequence of the liquid-like character of the 
post-cross-linking system on Icngthscales shorter than 
^L. As for the overall amplitude of H, this we expect 
to increase with (i) the fraction G of polymers that 
arc localized, (ii) the sharpness of localization, 1/£,l, 
(iii) the nematogen-nematogen aligning interaction J, 
and (iv) the length £ of the nematogens; and we ex- 
pect it to decrease with the "measurement temperature" 
T (see below for more on this concept), because ther- 



mal fluctuations tend to moderate any aligning forces. A 
complementary microscopic calculation [9| bears out this 
expectation, yielding Hk = Hq exp(— fc^f£/2), where the 
amplitude Hq = GV2(£/^i)Vr. 

In terms of these ingredients, we take as a model for 
the Landau-type free-energy cost F associated with the 
induction of local nematic order in the post-cross-linking 
system the form: 

F = i /" ((^i + /:fc2+Hk){QkQ-k} 

-2{(Yk + (r/r,)HkQi:)Q_k})- (3) 

Here, J-^ is shorthand for J d^ k / {2tt)^ , k^ is the squared 
length of the vector k, and the i?k is the Fourier trans- 
form / d^r R{r) exp(ik-r). In addition, curly brackets — 
as in {S S'} — indicate the trace of the product of the 
tensors S and S', i.e., J2dd'=i ^dd'S'^^,^. Furthermore, 
A characterizes the aligning tendencies of nematic free- 
doms; its value can be obtained by a complementary 
microscopic calculation as J'^/T (see Ref. Q), and C 
is the Frank constant for nematic order. Moreover, t 
is the reduced measurement temperature, defined to be 
(T — T*)/T* , where T is the measurement temperature 
(i.e., the temperature at which the system is maintained, 
in equilibrium, long after the cross-linking process), and 
T* is the spinodal temperature for the spatially homoge- 
neous isotropic-to-nematic transition. We also introduce 
the temperature Tp of the equilibrium state into which 
cross-links are instantaneously created, where p stands 
for preparation. The occurrence of two temperatures, T 
and Tp, stems from the fact that elastomers and related 
systems are characterized by not one but two statisti- 
cal ensembles. One, which we call the preparation en- 
semble , provides a statistical description of the random 
(non-equilibrating, unmeasured) freedoms Q" that char- 
acterize the local alignment immediately prior to cross- 
linking. The other describes the equilibrium state of the 
system long after cross- linking was done, via the statistics 
of the equilibrating variables Q; we call it the measure- 
ment ensemble. 

The free energy ^ consists of two terms. The first two 
elements of the first term constitute the familiar Landau- 
de Gennes free energy; they describe the free-energetic 
cost of inducing nematic alignment from the unaligned 
state. The second term incorporates what we have de- 
scribed above, viz., the influences of (i) the configuration 
of the rod-like constituents at the instant of cross-linking, 
via Q", together with (ii) the memory- independent ran- 
dom field Y caused by the localized polymers post cross- 
linking. The contribution to F/T involving Q'' carries 
a factor {J/Tj,){G(.^ /^lf{J/T). The two temperature 
factors show that the network is better able to store a 
given pattern Q*' when the preparation temperature Tp 
is lower and, similarly, better able to elicit Q" from Q the 
lower the measurement temperature T . Taken together. 



the two terms in F are minimized by the most probable 
nematic configuration Q, which is 



Qk = 



Yk + {T/Tp) Hk Q^ 
At + Ck'2 + Hk 



(4) 



The third element in the first term of the free energy ([3]) is 
a new and central element. It encodes the essential phys- 
ical difference between nematic fiuids and nematic elas- 
tomers, viz., the elastomer's possession of a network that 
is localized randomly and fluctuating thermally, and is, 
furthermore, liquid-like at sub-localization-length scales 
and solid-like at larger scales. As a result of this solid- 
ness at larger scales, creating a pattern of nematic align- 
ment (either by making relative translations of localized 
nematogens whilst maintaining their orientations, or by 
locally reorienting the nematogens) carries an additional 
free-energy cost, compared to nematic fluids, owing to 
the need to compete with localization forces. In a sim- 
ilar vein, by completing the square with respect to the 
first and second terms in Eq. ([3]), and thus arriving at 
the form 



{At + Ck" + Hk) { (Qk - Qk) (Q-k - Q-k) } (5) 



(up to a non-thermally fiuctuating term), the third cl- 
ement in Eq. ([5]) is seen to give a nonlocal free-energy 
cost of creating a departure from the nematic pattern 
Qk. This cost arises because the network mediates ad- 
ditional nematic-nematic interactions. Moreover, these 
nonlocal interactions are not independent of the memory- 
independent random field Y that acts to produce local 
nematic alignment, as this field also originates with the 
network. Thus, it is relatively costly, energetically, to 
create a nematic departure from Q that is essentially uni- 
form over a lengthscale rather larger than ^l, at which 
scale the solidness of the network becomes pronounced. 
On the other hand, the cost is relatively mild if the de- 
parture varies only over some lengthscale rather shorter 
than ^L , where the system has a more liquid- like charac- 
ter. 

Various averaged diagnostics of the system involving 
the local nematic order Q(r) may be considered via 
F. These averages come in two types. First, there are 
disorder-averaged quantities (denoted by [•••]), by which 
we mean quantities averaged over suitably distributed Y 
and Q*^. Second, there are thermally averaged quantities 
(denoted by (■••))' by which we mean quantities aver- 
aged over the measurement ensemble, i.e., the familiar 
statistical-mechanical ensemble of the equilibrating free- 
doms, in this case Q. Here, we focus on two particular 
diagnostics of nematic elastomers. The first, the thermal 
fluctuation correlator C^, defined via 



C^{r 



{(Q(r)-(Q(r)))(Q(r') 



(Q(r')))}) 
(6a) 



characterizes the strength of the thermal fluctuations of 
the nematic alignment away from the local mean value 
as well as the spatial range over which these fluctuations 
are correlated. Inter alia, through its range, C^ can be 
used to detect the occurrence of a continuous phase tran- 
sition. The second is the glassy fluctuation correlator C^, 
defined via 



C«(r-r')^[{(Q(r))(Q(r'))} 



(6b) 



This is a diagnostic of any randomly frozen (i.e., time- 
persistent) nematic order present. In particular, its value 
at the origin, C'~^{r)\r=o, is the nematic analog of the 
Edwards- Anderson order parameter for spin glasses [lO[ , 
and measures the magnitude of the local frozen nematic 
order; hence the name glassy correlator. Moreover, how 
C^{r) varies with the separation r determines the spatial 
extent of regions that share a roughly common nematic 
alignment. For isotropic liquids the glassy correlator au- 
tomatically vanishes, but for IGNEs it is nonzero. The 
correlators C^ and C*^, which diagnose the measurement 
ensemble, give valuable information about the physical 
structure of nematic elastomers [11|. (Other correlators, 
including ones that, in addition, probe the preparation 
ensemble, can also be considered. A broad account of the 
issues that result from the presence of a multiplicity of 
ensembles and the correlators that diagnose them will be 
given in a forthcoming paper J12l|.) 

One could also consider the disorder-averaged quan- 
tity [(Q(r"))]; it, however, vanishes, owing to the macro- 
scopic isotropy of the post-cross-linking state. On the 
other hand, the thermal average of the local order pa- 
rameter for a specific realization of the quenched disor- 
der (Q(r)) is maintained at a nonzero, time-persistent, 
random value, which we shall compute shortly. This 
nonzero value is the result of the partial trapping, by 
the network, of the orientational randomness Q*' present 
at the instant of cross-linking, together with the memory- 
independent random field Y of the network, post cross- 
linking. The free energy ([3]) is quadratic in Q, and there- 
fore the computation of (Q) and C^ using the statistical 
weight exp(— F/r) is elementary, yielding 



(Q-k) = 



Yk+(r/rp)HkQ^ 



At+Ck^+Ui^ 

Qk-(Qk))(Qk'-(Qk' 



TfiDSk+k'.O 

At + £fc2 + Hi, 



(7) 
.(8) 



Here, fio = (-D — !)(£* + 2)/2 counts the number of de- 
grees of freedom of Q and takes the value 5 for D = 3. 
Note that we have chosen units in which Boltzmann's 
constant has the value unity. 

To perform the average over the quenched random 
variables Y and Q° we must adopt a model for their 
statistics that is consistent with the physical origin each 
has. The choice we make is that Y and Q" arc indepen- 
dent, Gaussian-distributed random fields, with vanish- 



ing means and non-vanishing variances, the latter being 
given by 



[{Q°kQ"k'}] 

[{YkYk.}] 



T 



AH^ + C^k^ ' 



rHk(5k+k',o 



(9a) 
(9b) 



Here, A'^ is the preparation-ensemble counterpart to A, 
and a complementary microscopic calculation shows that 
its value is given by J'^/Tp (see Ref. [9|]). Similarly, C^ is 
the preparation-ensemble counterpart to C. The statis- 
tics of Q° depends on tp [i.e., the parameter that en- 
codes the temperature Tp of the preparation ensemble 
via tp = (Tp — T*)/T*]; it does not depend on H, be- 
cause H encodes the physics of random but imperfect 
spatial localization, and this only comes into being post 
cross-linking. (The impact of Q*^ does depend on H, as 
H controls the relaxation of Q from Q" to its equilibrium 
value, post cross- linking.) By contrast, the statistics of 
Y does depend on H; this is because H characterizes the 
typical value of the memory-independent random field 
that results from the random imperfect spatial localiza- 
tion of the polymers constituting the network. In view 
of their distinct origins it is natural that Y and Q° are 
statistically uncorrelated. However, it is also natural, 
at least hcuristically, that the H that characterizes the 
"orientational caging" induced by the network (via Y) is 
the same H that determines the fidelity with which the 
network preserves the orientational order present imme- 
diately post cross-hnking (i.e., Q°). It is natural because 
localization that is sharper and more widespread (i.e., 
involves a larger localized fraction) both creates more in- 
tense network-induced orientational caging and enhances 
the trapping-in of the local nematic order present im- 
mediately post cross-linking. Our physical expectation, 
borne out by a complementary microscopic analysis (see 
Ref. [9(]), is that such strengthening of the localization 
would more strongly enhance memorization than orien- 
tational caging. This expectation is consistent with the 
phenomenological choice presented here, in which the 
corresponding contributions to the random anisotropy 
field, Eq. ^, scale as -s/H for the caging (i.e., Y) part 
and H for the "memorization" (i.e., Q°) part. 

Returning to the disorder-averaged diagnostics — the 
mean value [(Q)] and the correlators C^ and C'^ — we 
complete their computation by using the statistics of the 
quenched disorder, Eqs. (jH]), to arrive at 



[(Qk)]=0, 
Ck = TfiD 



At + Ck^ + Hk ' 



Ct = TfiD ■ 



{At+Ck^+Ri^f ■ 



(10a) 
(10b) 

(10c) 



Having computed the correlators C'^ and C^, we now 
set about using them to study how the presence of a 



network modifies the organizational behavior of nematic 
freedoms. To do this, we first note that there are two 
emergent lengthscales present in IGNEs: (i) the typical 
localization length, ^^ , quantifying the sharpness of local- 
ization of polymers belonging to the network, and (ii) the 
intrinsic nematic correlation length, ^n [= \/^/M\, de- 
scribing the range over which nematic freedoms would be 
correlated if there were no network present. On the other 
hand, we have the strength of the memory-independent 
random field Y, which is characterized by VHo. In what 
follows, we shall study the effects of these three parame- 
ters on C'^ and C^ at two preparation temperatures, one 
far above T* and the other just above it. 



Disorder strength 



t,G,o 



Weak (Ho < H, 



2 l-(Ho/Hc) 
N l+{Ho/^t) 



CXJ 

2?i ~r ^?Af l+{Ho/>lt) 



THo/ 



SAtZJ'iL 



Strong (Ho > H^ 



iCi/ln(Ho/Hc) 

~ci/(i+4) 

iCi/ln(Ho/H,) 

~ 2« 



Tf^IWZ^ 



TABLE I: Values of the correlation lengthscales {^T,d and 
^G,d), the oscillation wavelengths {^t,o and ^g,o), and the 
intensity of local nematic alignments [C (r = 0)] in the weak- 
and strong-disorder regimes for the case of IGNEs crosslinked 
at temperatures far above T*. 



First wc consider the behaviors of C and C for sys- 
tems prepared at some temperature far above T* , so that 
any local nematic order present immediately post cross- 
linking (and thus available for trapping in) is spatially 
correlated only over distances far shorter than the typi- 
cal localization length ^l ', see Table [J This implies that 
the local nematic order arising from Q° would be heav- 
ily "washed out" by thermal fluctuations of the network. 
Thus, in this regime, the dominant contribution to the 
trapped-in local nematic order originates in the memory- 
independent random field, Y, and concomitantly, the 
contribution arising from Q° may be neglected. 

C^ and C'^ exhibit qualitatively distinct behaviors in 
two regimes, depending on the strength of the random 
field. For Hq < He (where He = 2£/^£ — the weak dis- 
order regime), C'^ and C simply decay with increasing 
real-space separation. More specifically, by examining 
their small wave-vector behaviors we ascertain that the 
respective associated correlation lengths S,T,d and ^G,d 
have the values given in Table [J We see from the be- 
havior of £,T,d the physically reasonable result that the 
random network, with its thermal fluctuations, serves 
to shorten the nematic thermal fluctuation correlation 
length from the value it would have absent the network. 



a phenomenon that a conventional (i.e., non-thermally 
fluctuating) random-field approach would not capture. 
As for ^^ ^, it comprises two parts. One (oc Crd) arises 
from the nematic thermal correlations; the other (oc ^£) 
comes from the local aligning effect exerted by the cage. 
The fact that £,G,d increases with S,l does not mean that 
a more weakly cross-linked network (for which ^l would 
be larger) aligns the ncmatogens more effectively. Whilst 
the lengthscale of aligned regions ^G,d may increase, the 
magnitude oiC'^ , which governs the intensity of the align- 
ment locally, decreases, as can be seen from Table HI 

On the other hand, the simple decay of the correlators 
at weak disorder can give way to oscillatory decay at 
strong disorder, i.e., Hq/Hc > 1. C^ has such behavior 
regardless of T, whereas C'^ has it only for sufficiently 
small T. The oscillation wavelengths £,t/g.o are deter- 
mined via the radii of the shells in wave-vector space on 
which the corresponding correlators are maximal. Thus, 
we arrive at an explicit (and, notably, T-independcnt) 
formula S,t,o = ^l/-\/2 ln(Ho/Hc) and an implicit one for 
^G.o, viz.. 



1 + {^N/^G,of + ^i^N/^Lf - (HoMi)e 



(11) 

The value of S.T/G.d in this strong-disorder regime, given 
in Table HI is estimated via the widths of the peaks of 
C'^l'^ . Upon decreasing ^l at fixed ^n , the value of i,T,d 
tends to ^^ from above, indicating that the network is 
limiting the range over which the thermal nematic fluc- 
tuations are correlated; on the other hand, i,G,d remains 
at the scale of ^^ , indicating that the range of coherent 
nematic alignment is circumscribed by the network's typ- 
ical localization length. Concomitantly, there is a growth 
in the intensity of the local alignment C'^(v = 0). 

Having considered the behaviors of C"^ and C'-^ for sys- 
tems prepared at high temperatures, we now consider the 
corresponding behavior for systems prepared at a tem- 
perature just above T* , so that the local nematic order 
present immediately post cross-linking is spatially corre- 
lated over distances larger than the typical localization 
length of the elastomeric network. As one can see from 
Eq. (jlObp . the behavior of C"'" is unchanged, undergoing 
simple decay in real space at weak disorder and oscil- 
latory decay at strong disorder. Conversely, C^ exhibits 
behavior qualitatively different from that of a system pre- 
pared at high Tp, because it now receives its dominant 
contribution from the memorization of Q°. Specializing 
tot 
given by 



tp, we sec from Eq. (jlOcp that C^ is approximately 



c2 



A^o TTT 



T 



T 



TpJ A% + COk^' 



(12) 



which is proportional to the correlator of the thermal ne- 
matic fluctuations immediately post cross- linking. This 
indicates that the pattern of these thermal fluctuations 
has been faithfully memorized by the network. 



To conclude, we have proposed a physically motivated 
phenomenological model to describe isotropic-genesis nc- 
matic elastomers. Our approach enables us to charac- 
terize the random nematic alignment present, post cross- 
linking, in terms of both the partial memorization by the 
network of the nematic alignment present at the genesis 
of the network and the structural anisotropy presented 
by the network itself. Our approach also reveals that 
networks that have sufficiently strong aligning capabili- 
ties exhibit a novel, oscillatory-decaying type of spatial 
correlation of nematic alignments. As we shall explore 
in detail in two companion papers [9|, ll6| , the physical 
origin of these oscillatory correlations is the thermally 
compliant nature of the random fields in these systems. 
The nematogens subject to similar random anisotropic 
environments tend to aggregate whilst those subject to 
different anisotropic environments tend to segregate. In 
addition to its specific results, the present work un- 
derscores the necessity of broadening the concept of a 
quenched random field to incorporate not only the tra- 
ditional, "frozen" type, which does not fiuctuate ther- 
mally [J, [sl, Il3l4l5 1 , but also the type in which the frozen 
character of the random field is present only at longer 
lengthscales, fading out as the lengthscale progresses to 
below a characteristic localization length, owing to the 
thermal position fluctuations of the network constituents. 
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of the correlators, it is straightforward to reconstruct the 

full structure of the corresponding fourth-rank Cartesian 

tensors, by appending suitable isotropic tensor factors 

constructed from Kronecker deltas. 
[12] X. Xing, B. S. Lu, F. Ye and P. M. Goldbart, Manuscript 

in preparation (2012). 
[13] Y. Imry and S. K. Ma, Phys. Rev. Lett. 35, 1399 (1975). 
[14] L. Radzihovsky and J. Toner, Phys. Rev. B 60, 206 

(1999). 
[15] D. E. Feldman, Phys. Rev. Lett. 84, 4886 (2000). 
[16] F. Ye, B. S. Lu, X. Xing and P. M. Goldbart, Manuscript 

in preparation (2012). 



[1] M. Warner and E. M. Terentjev, Liquid Crystal Elas- 
tomers. Oxford University Press (2003). 

[2] K. Urayama, Macromol. 40, 2277 (2007). 

[3] K. Urayama, E. Kohmon, M. Kojima, and T. Takigawa, 
Macromol. 42, 4084 (2009). 

[4] S. V. Fridrikh and E. M. Terentjev, Phys. Rev. E 60, 
1847 (1999). 

[5] L. Petridis and E. Terentjev, Phys. Rev. E 74, 051707 
(2006). 

[6] N. Uchida, Phys. Rev. E 62, 5119 (2000). 

[7] J. V. Selinger and B. R. Ratna, Phys. Rev. E 70, 041707 
(2004). 

[8] P. M. Goldbart, H. E. Castillo and A. Zippelius, Adv. 
Phys. 45, 393 (1996). 

[9] B. S. Lu, F. Ye, X. Xing and P. M. Goldbart, Manuscript 
in preparation (2012). The microscopic model exam- 
ined in this paper involves a system of dimers that 
are randomly and permanently connected by harmonic 
springs, and interact with one another via Maier-Saupe 
and excluded- volume forces. 
[10] S. F. Edwards and P. W. Anderson, J. Phys. F: Metal 

Phys. 5, 965 (1975). 
[11] Although we are focusing on the Cartesian scalar aspects 



